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We investigate the binding energy per nucleon and hadron properties in infinite and homogeneous 
nuclear matter within the framework of the in- medium modified Skyrme model. We first consider the 
medium modifications of the single hadron properties by introducing the optical potential for pion 
fields into the original Lagrangian of the Skyrme model. The parameters of the optical potential are 
well fitted to the low-energy phenomenology of pion-nucleus scattering. Furthermore, the Skyrme 
term is also modified in such a way that the model reproduces bulk properties of nuclear matter, 
in particular, the binding energy per nucleon. The present approach is self-consistent: the single 
hadron properties in a nuclear medium, their effective in-medium interactions, and the bulk matter 
properties are treated on the same footing. 
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I. INTRODUCTION 



The equation of state (EOS), which gives the density dependence of the binding energy per nucleon for a given 
nucleus, has been one of the most importance issues in nuclear many-body problems. There is a great amount 
of different theoretical approaches in trying to describe the EOS and, clearly, we can mention only some of few 
representatives [3-0] ■ In general, those approaches and corresponding representatives can be classified into three 
classes as microscopic many-body approaches , effective field theories [H, HJ , and phenomenological methods 0] ■ 
These works for many-body problems provide very useful tools for understanding properties of dense and hot matter. 

On the other hand, the Skyrme model 0, Q also presents a simple but good framework for investigating bulk 
properties of nuclear matter. In general, one can classify various Skyrmion approaches into two subclasses: In the 
first one investigations are mainly devoted to the classical crystalline structure and its behavior under the extreme 
conditions [HI El • In the second approach the properties of exotic many-baryon systems were treated (l2l - fl4| . There 
are also some early attempts to explain many-body systems within the Skyrme model (l5l [HJ , considering the single 
skyrmion in hypersphere. 

Moreover, there is an another alternative way to examine properties of nuclear matter within the Skyrme model, 
i.e. to study the properties of the single skyrmion in nuclear matter [l7| . Furthermore, the medium- modified Skyrme 
model [13, El , in connection with quantum-mechanical variational methods [l9[ (like the Hartree-Fock method) , can 
be applied to the analysis of the bulk properties of nuclear matter [2(j. To perform this analysis, it is essential to know 
the properties of the single hadron and the NN interaction in a symmetric [13, Gil as weu as asymmetric [2l|) nuclear 
environment. Consequently, the behavior of the hadrons in nuclear medium must be taken into account. However, 
the variational calculations can only estimate the upper boundary value of minimized quantities. Thus, one may still 
consider how to improve the results. 

As a more realistic approach, the in-medium modified Skyrme model [l7| itself can be used so as to reproduce 
the properties of the single hadron in nuclear matter as well as those of matter in bulk. This is the aim of the 
present work. To carry out the proposed goal, we consider not only the changes of the kinetic and mass terms of the 
standard Skyrme Lagrangian as done in Ref. [l7| but also possible modifications of the Skyrme term. It is well known 
that Skyrme's quartic stabilizing term may be related to vector mesons [22j that can be realized in implicit gauge 
symmetry of the nonlinear sigma model Lagrangian |23j . In this sense, the modification of the Skyrme parameter 
may be pertinent to the changes of the vector mesons in nuclear matter. 

The present work is organized as follows: in the next Section, we explain how the original Skyrme model can be 
modified in medium. In Section III, we present the corresponding numerical results and discuss them. In particular, 
we show that the quartic Skyrme term can be modified in such a way that it minimizes the binding energy per 
nucleon. We discuss the changes of the mass splitting of the nucleon and the A isobar and those of the irNN coupling 
constant. We also estimate the symmetry energy in the Weizsacker-Bethe-Bacher formula. The last Section is devoted 
to summary and outlook of this work. 



II. MEDIUM MODIFICATION OF HADRON PROPERTIES 

In Ref. [I?}, the in-medium modified Skyrme Lagrangian was presented, the mass term being modified based on 
the phenomenology of low-energy pion-nucleus scattering. The modified mass term leads also to the changes of the 
kinetic term. In the present work, we additionally consider the modification of the Skyrme stabilizing term. The 
resulting Lagrangian is given as follows: 

F 2 f f)TI\ /8TI^\ F' 2 1 

c = tH*) (Tr)-it°»M Tr < w >-< w, > + 3^ 1Y i t,,i '» c '' c,, * t, i 2 

F 2 m 2 

+ ^^a s (r)Tr(C/ + [/t _ 2 ) , (1) 
16 

where F^ denotes the pion decay constant, e is the Skyrme parameter, and m w stands for the pion mass. The medium 
functionals a s and a p are written in the following forms 



4irb p{r)f 4irc p(r) 
a s = l j , a p = l- — — — rT , (2) 



which represent the influence of the surrounding environment on the properties of the single skyrmion. These param- 
eters are related respectively to the corresponding phenomenological S- and P-wavc pion-nucleus scattering lengths 
and volumes, i.e. bo and cq, and describe the pion physics in a baryon-rich environment [24j . The density of the sur- 
rounding nuclear environment is given by p, g' denotes the Lorentz-Lorenz or correlation parameter, / = l + m^/m 1 ^ 00 
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represents the kinematical factor, and m^° e is the nucleon mass in free space. In addition to these changes of the 
mass and kinetic terms done in Ref. fl7| . we introduce the new density-dependent functional 7(1*) = j(p(r)) which 
provides the in-medium dependence of the Skyrme parameter, i.e. e 2 — > e* 2 = e 2 7. As mentioned in Introduction, 
this medium modification can be related to the density dependence of vector meson properties. Furthermore, to fix 
this additional dependence in the present work, we will concentrate on the bulk properties of infinite and homogenous 
nuclear matter with constant density (p = const). Note that the Lagrangian {T]) is modified in such a way [l7T | that at 
zero density it reduces to the original Lagrangian of the Skyrme model and at the linear approximation it reproduces 



the well-known equation of the pion fields in nuclear medium [24 1 . 

In order to analyze the field equation for the classical pion field in homogenous nuclear matter one can choose the 
spherically symmetric "hedgehog" form for the boson field U = exp{m • rF(r)}, where n denotes the unit vector 
in coordinate space, r are the usual Pauli matrices, and F(r) stands for the profile function of the pion field. The 
pertinent field equation, which is given as 



F"(x) (a pl x 2 + 8s 2 ) + 2a pl xF'(x) - a s -/P 2 x 2 s + (^F'(x) 2 - a p 7 - sin(2F) = 0, 

is obtained by minimizing the medium-modified mass of the static skyrmion 



(3) 
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In the last two expressions, we have introduced the dimensionless variable x = eF^r and the new constant ft = 
m w /(eF 7r ). 

The collective quantization of the classical skyrmion |9j yields the in-medium modified nucleon mass and the 
corresponding A — N mass splitting respectively as 

3 3 
m* N = M£ + —, m* A _ N = — , 

where A* denotes the in-medium moment of inertia of the skyrmion. The meson-baryon vertices in nuclear matter 
can be derived by calculating the in-medium modified irNN form factor [l8| 

2 . _ 4ttM* f°°h(qx) „ , , 3 



Gl NN {q ) = „ 2ff a p / Sk(x)x Ax, (6) 

•JC -T-K JO 

where q = q/eF^, ji(qx) is the spherical Bessel function with order 1 and S n (x) is defined as 

S w {x) = -(2x- 1 F' + F") cos F+(F' 2 + 2x-' 2 + a s a p 1 ml) sin F. (7) 

Using the Lagrangian given in Eq. ([T]) , one can calculate the in- medium modifications of the single nucleon properties 
and the pion-nucleon coupling constant. Note that the parameters of the model are fitted to be F^ = 108.78 MeV and 
e = 4.85 so as to reproduce the experimental values of the nucleon and A in free space. Consequently, the pion mass 
is also fixed to be its experimental value for the neutral pion, i.e. = 134.98 MeV. A set of values of parameters in 
the medium functionals ([2]) are taken from the analysis of phenomenological data for pion-nucleus scattering (24[. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



As is clear from the discussions in the previous Section, if the modification of the Skyrme term is ignored, the values 
of all input parameters arc fitted to the phenomenology or taken from it. When, however, the Skyrme term is modified, 
we introduce one additional functional 7(p) which, in general, may be related to the vector meson properties in nuclear 
matter. One can also note that the lessening value of the Skyrme parameter in nuclear medium may correspond to 
a decrease of the gp^T coupling and, therefore, to the change of the rho meson width in nuclear matter or to a 
diminishing value of its mass in medium, i.e. m*/m p < 1. There are experimental indications to those changes of the 
p meson properties [25U27I ] and the theoretical predictions [H, [29| . Following the ideas presented in those theoretical 
approaches, one may be able to fit 7. However, it is still under debate how the properties of the p meson undergo in 
medium both theoretically and experimentally, and mo del- dependent. Thus, in the present work, we will proceed to 
fit the form of the functional 7 to the bulk properties of nuclear matter rather than following a specific model. 
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As a first step, wc explicitly choose its form to reproduce the first coefficient (volume term) in the semiempirical 
Weizsackcr-Bethc-Bacher mass formula. Then the binding energy per nucleon at a given density can be defined simply 
as 

AE B=1 = m* N (p) - m% cc . (8) 

This is somehow a crude approximation but a comprehensive one. One could even fit the form of 7 by investigating 
different terms in the mass formula and by examining the interplay between them. However, within the present work, 
the approximation defined in Eq. (JSJ will be enough for the qualitative analysis of the changes due to the modification 
of the Skyrme parameter. 

In the present calculation we have tried various forms of the dependence of 7 on the nuclear density p such as linear, 
quadratic, polynomial, exponential forms, etc. It turns out that the best fit to the ground state of nuclear matter is 
achieved by the following form 

7 (p)=expf--^0 , (9) 



1 + 7denP, 

where 7 num and 7den are variational parameters. Let us first discuss the behavior of the binding energy when the 
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FIG. 1. (Color online). The binding energy per nucleon as a function of p/po- The solid curve (left scale) corresponds to the 
parametrization of 7 in Eq. with 7 nU m = 2.1?n~ 3 , 7a en = 1.45m~ 3 and P-wave scattering volume Co = 0.21m~ 3 used. The 
dashed one (left scale) draws the case when 7 nU m = 0.8m^ 3 , 7d on = 0.5m^ 3 and P-wave scattering volume Co = 0.09m^ 3 . 
5-wave scattering length is fixed at bo = — 0.024m^ 1 and the correlation parameter has value g' = 0.7. The dotted one (right 
scale) shows the case that Skyrme term intact in nuclear matter and consequently j(p) = 1. The normal nuclear density is 
given as po = 0.5m 3 . 



Skyrme term is intact in nuclear matter, i.e. j(p) = 1. The corresponding binding energy is depicted as the dotted 
curve in Fig. [T] with the energy scale drawn at the right vertical axis. One can note that in this case the binding 
energy monotonically falls off as the density increases. In the language of the single skyrmion, it indicates that the 
skyrmion swells to a larger volume and all skyrmions of the system start to overlap. Thus, the density of the system 
continuously increases. This is not surprising, because the medium modification in this case can be simply related to 
that of the pion decay constant F n — > F* = F^^Ja^. For the moment, one can ignore the explicit chiral symmetry 
breaking term in the Lagrangian, because its influence to the stability is rather small in comparison with the effects 
coming from the first two terms. The decreasing value of the pion decay constant changes the contribution from the 
nonlinear kinetic term. As a result the skyrmions swell to the larger volume and it is necessary to prevent this by some 
mechanism. It implies that one must introduce either strong repulsive NN interactions at short distances or some 
mean-field mechanism as in the Walecka model @, 0). However, one interesting way to avoid this collapse may be to 
modify the Skyrme term and it is also physically motivated as discussed at the beginning of this Section. Moreover, 
it is much simpler and transparent to consider the modification of the Skyrme term. 
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We also present two different results with the modified Skyrme term in Fig. [TJ The solid and dashed curves draw 
the parametrization of Eq. © with the energy scale depicted at the left vertical axis. Note that the values of the 
variational parameters, 7 nU m and 7den, are chosen in a such way that the minimum of the binding energy occurs at the 
normal nuclear matter density and reproduces correctly the first coefficient in the empirical formula for the binding 
energy. The difference between these two curves is due to the fact that we use two different values of the P-wave 
scattering length, i.e. Co = 0.21m~ 3 for the solid curve and Co = 0.09m~ 3 for the dashed one. The results show that 
the dependence on the density is rather insensitive to the changes of input parameters from pion-nucleus scattering 
phenomenology and moreover the effect of the changes in 60 is even milder. 

In order to see the validity of Eq. ([9]), it is of great importance to examine the changes of other physical observablcs. 
Let us first discuss thermodynamic properties of nuclear matter. The pressure is given by the following formula 



V = P 



de 
dp 



P 



,dAEf 



dp 



(10) 



where e is the total binding energy of nuclear matter per unit volume. It vanished naturally at the equilibrium point 
and for the parametrization of Eq. (|9]). We want to emphasize that the pressure is always decreasing with the Skyrme 
term intacted. 

Another important quantity is the compressibility of nuclear matter expressed as 



A" 



Wo 
dp 2 



2 d 2 AE B=1 



dp 2 
1 

1 



9p§ 



P=P0 



d 2 a s 



ttP 



dp 2 
dx ( F 



e 

sin z F 



dx 



p=po ~ 
2 

a 



x 2 F' : 



sin F 



(11) 



2x 2 



sin 2 F - 



dx : 



F 



sin 2 F 



sin 2 F 



and the corresponding results are listed in Table U with two different values of the S- wave scattering length 60 [HI 
used. The results show that the compressibility of nuclear matter and the effective A — N mass difference are rather 



60 [m^ x ] 


co ["O 


7num [m„ 3 ] 


7don [77V 3 ] 


K [MeV] 


m* N _ A [MeV] 


-0.024 


0.21 


2.098 


1.451 


1647.47 


105.21 


-0.024 


0.15 


1.448 


0.998 


1148.18 


129.39 


-0.024 


0.09 


0.797 


0.496 


582.79 


170.34 


-0.029 


0.21 


2.106 


1.506 


1637.16 


107.13 


-0.029 


0.15 


1.444 


1.031 


1142.00 


131.59 


-0.029 


0.09 


0.785 


0.502 


580.03 


172.91 



TABLE I. Compressibility of nuclear matter K and an effective A-nucleon mass difference m* N _ A at the normal nuclear matter 
density po = 0.5m 3 . The variational parameters 7 nU m and 7d cn are fitted to reproduce the minimum of the binding energy per 
nucleon ~ 15.7 MeV at the normal nuclear matter density. The correlation parameter is taken to be g' = 0.7. 

stable under the change of bo. On the contrary, they are quite sensitive to the value of P-wave scattering volume Cq. 
At the empirical value of cq = 0.21m~ 3 , the compressibility turns out to be very large (A ~ 1640 MeV) in comparison 
with those obtained in rclativistic Dirac-Brueckner-Hartree-Fock approaches 0, H| and in the Walecka model p} • We 
find that as lower values of Co are used K is noticeably decreased. For example, for c = 0.09m~ 3 the compressibility 
is already consistent with that of the Walecka model (A' ~ 580 MeV). If one uses even a smaller value of Co such 
as Co = 0.06m~ 3 , the result of A is further brought down to be comparable with that in Dirac-Brueckner-Hartree- 
Fock approaches (A ~ 300 MeV), which is close to the empirical value. It indicates that the present work prefers 
smaller values of cq than that used in the pionic atom analysis as far as the compressibility is concerned. Note that 
a similar conclusion about Co was drawn from the analysis of an effective axial-vector coupling within the original 
medium-modified Skyrme model (l7j . We remind that K is sensitive to the position of the saturation point. Fitting 
the saturation point at slightly lower densities, we see that the compressibility decreases drastically. However, the 
situation may change if one considers a more accurate approximation with the surface and symmetry energy terms 
explicitly taken into account in Eq. ([5]). 

In Fig. the dependence of the A — N mass difference on the nuclear matter density is drawn. The results show 
that the modified Skyrme term leads to rather different results from that without the modification. With the Skyrme 
term modified (see solid and dashed curves), the results of m^_ N fall off faster as the density increases in comparison 
with that with the original Skyrme term (see dotted curve). Of course, this is due to the explicit density dependence 
of the moment of inertia ([5]) through the additional density functional 7. It implies that it is easier to make the 
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300 




FIG. 2. (Color online). The density dependence of the A — N mass difference in nuclear matter. The notations and input 
parameters are the same as those in Fig. [T] 

nucleoli excited to the A state in nuclear matter, which seems more realistic than that without the modification of 
the Skyrme term. 

In Fig. [31 the changes of the irNN coupling constant are depicted. Here the modifications of the Skyrme term bring 




FIG. 3. (Color online). The dependence of ttNN coupling constant on the density. Notations are similar to those in Fig. [T] 



about more dramatic results. When the Skyrme term is intact, i.e. 7 = 1, g^NN monotonically decreases as the density 
increases (see the dotted curve in Fig. [3J. However, when one introduces the density dependence of the Skyrme term, 
the results are noticeably changed. Using the paramctrization of Eq. ([9]) with the different values of input parameter 
Co, we find that the ttNN coupling in medium changes drastically. For example, with the value of Co = 0.09m~ 3 the 
in-medium pion-nucleon coupling constant g* NN starts to increase monotonically up to high (p ~ 5po) densities as 
drawn in the dashed curve. The g^Ar with cq = 0.09m~ 3 will disappear at around p ~ 5po- On the other hand, if one 
uses Co = 0.21m~ 3 , it is getting increased up to the normal nuclear matter density and stays more or less constant. 
Then it slowly falls off as the density increases. When it approaches the critical point p ~ 2.54po, it drops sharply 
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and goes to zero (see the solid curve in Fig. [3]). Above the critical point (p > p a i t ~ 2.54po), the skyrmion does not 
exist. It is not surprising, because the sign of the combination a p 7 in Eq. is changed and therefore there is no 
stable solitonic solution anymore. 

Let us draw the attention again to the bulk properties of nuclear matter in order to understand the modifications 
of the Skyrme term better. Following Klebanov [lOj . after quantizing and using the formula presented in Ref. [30| . 
one can estimate the symmetry energy in the semiempirical formula for the nuclear binding energy: 

-Bsym = ^ m*A-N ■ ( 12 ) 

This crude formula of the symmetry energy already provides the enlightening results. For example, 7 parameterized 
as in Eq. ©, E sym (p ) ~ 14.19 MeV for c = 0.09m~ 3 whereas E sym (p ) « 8.71 MeV for c = 0.21m,,: 3 . These 
results for the symmetry energy must be compared with the experimental one E sym ~ 20 — 30 MeV. The order of 
the symmetry energy calculated within the Skyrme model is comparable to the experimental data. To estimate the 
symmetry energy more accurately, however, one should consider the minimization of the whole binding energy taking 
into account the interplay between the different terms in the mass formula as we stated already. Moreover, one should 
consider the effects of finite nuclei and explicit isospin-breaking effects and so on. We want to mention that this is also 
possible within the in-medium modified Skyrme model and can be done as in Refs. (2ll. l3ll] . an additional modification 
of the Skyrme term being performed as was done in the present work. 



IV. SUMMARY AND OUTLOOK 



In the present investigation, we aimed at studying the modifications of the quartic term in the Skyrme model. The 
results from this work shows that the additional modifications change dramatically the whole picture and allows one 
to understand the role of the modifications in a more comprehensive way. One can note that an alternative approach 
to many-body systems within the Skyrme model fl2l [l3| points to the changes of the input parameters (so called 
"calibration" ) according to the number of baryons in the system. Within our approach these changes were shown in 
a more realistic and transparent way and were treated not only at the level of the system but also at the level of its 
constituents. 

From the previous studies we know that the large renormalization of the nucleon mass in nuclear medium causes 
one of the difficulties to bind the infinite nuclear matter [20[ and to reproduce the correct values of the Nolen-Schiffcr 
anomaly in mirror nuclei within the in- medium modified Skyrme model plj . The relatively small change of the 
nucleon mass in nuclear matter within the present approach is an interesting result, which can be used to reproduce 
the correct value of the Nolen-Schiffer anomaly. In addition to the small nucleon mass renormalization, the dramatic 
changes in the pion-nuclcon coupling constant gives the opportunity to revise the previous investigation on nuclear 
matter related to the quantum- mechanical many-body problems [201 ] . 

More completely, one can improve the present approach by treating the density of the system in a fully consistent 
way. Since it is a sum of the single skyrmion densities in a given initial configuration and one can study the equilibration 



of the matter to its ground state, taking into account the deformation effects on the properties of its constituents [31 



Having considered the modification of the Skyrme term, we assert that the present approach is self-consistent: 
it treats the single hadron properties, the hadron-hadron interactions and the bulk matter properties on an equal 
footing. Moreover, it is closely related to phenomenological low-energy data at the single-hadron level as well as at 
the level of hadronic systems. The present approach can be extended to the studies of the properties of finite nuclei 
and its constituents. 
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